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ABSTRACT 

A  model  forecasting  the  wave  height  and  the  longshore  current  distribution  inside 
the  surf  zone  based  on  the  formulations  by  Thornton  and  Guza  (1983),  (1986),  are  ap- 
plied to  an  extensive  set  of  both  laboratory  and  field  data  for  the  purposes  of  testing  and 
modification.  The  models  are  tested  on  planar  beaches  as  well  as  barred  beaches  for  a 
variety  of  wave  conditions. 

The  wave  transformation  model  is  based  on  solving  the  energy  flux  equation  using 
a  bore  dissipation  mechanism  and  describing  the  random  wave  heights  with  the  Rayleigh 
distribution.  The  two  model  parameters  B  and  y,  where  B  describes  the  amount  of  foam 
of  a  breaking  wave  and  y  is  the  proportionality  constant  which  relates  the  rms  wave 
height  to  the  water  depth,  are  combined  into  a  single  parameter  BG.  The  combined 
parameter  BG  is  shown  to  be  a  function  of  deep  water  surf  similarity  parameter.  Ap- 
plied to  the  present  data  sets,  the  rms  error  of  the  measured  wave  height  and  the  model 
predicted  wave  height  was  usually  less  than  9°/ b  and  ranged  from  1.5%  to  15.7%  with 
a  mean  of  5.3%  and  a  standard  deviation  of  3.1%  for  the  whole  74  data  sets.  The 
wave  transformation  model  is  highly  robust  in  describing  the  wave  height  distribution  in 
the  surf  zone. 

The  longshore  current  model  is  based  on  solving  the  steady  state,  alongshore  mo- 
mentum balance  for  straight  and  parallel  contours  using  the  radiation  stress  concept. 
The  model  requires  specifying  the  bed  shear  stress  coefficient  (cf)  and  turbulent  mixing 
coefficient  ( N  >.  Applied  to  the  present  data  sets  the  rms  error  between  the  measured  and 
modeled  longshore  current  values  ranged  from  4.5%  to  55.5%  with  a  mean  of  24. 6%. 
The  turbulent  mixing  is  not  required  for  planar  beaches,  but  it  is  required  for  barred 
beaches  to  describe  the  longshore  currents  inside  the  surf  zone.  The  mean  value  of  the 
bed  shear  stress  coefficient  for  the  present  longshore  current  data  is  0.006  and  the  mean 
lateral  mixing  coefficient  is  0.006. 
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I.     INTRODUCTION 

A.     BACKGROUND 

The  wave  height  and  longshore  current  distribution  inside  the  surf  zone  is  a  primary 
concern  of  Coastal  Engineering  and  Naval  Operations.  These  distributions  are  impor- 
tant for  amphibious  operations  both  at  the  decision  making  stage  and  during  the  beach 
assault  in  order  to  reduce  the  percentage  of  casualties.  The  type  of  landing  craft  to  be 
chosen  depends  on  the  wave  and  current  conditions  in  the  surf  zone. 

As  waves  move  from  deep  water  to  shallow  water,  they  are  transformed  due  to  the 
effects  of  shoaling,  refraction  and  dissipation.  Shoaling  results  in  decrease  in  wave 
length  and  an  increase  in  wave  height.  Wave  refraction  results  in  change  in  wave  direc- 
tion. Shoaling  and  refraction  are  reasonably  well  defined  by  theory.  Dissipation  is  due 
to  bottom  friction  and  wave  breaking.  Wave  breaking  is  the  dominant  dissipative 
mechanism  inside  the  surf  zone.  At  a  certain  depth  the  waves  become  unstable  and 
break  in  either  spilling,  plunging  or  surging  form,  depending  on  the  deep  water  wave 
steepness  and  the  beach  slope.  After  breaking  the  wave  height  decreases  as  the  energy 
is  converted  into  turbulence  and  finally  lost  to  heat.  If  the  waves  approach  the  beach 
at  an  angle  then  they  cause  a  longshore  current  parallel  to  the  shoreline  within  the  surf 
zone  because  of  the  change  in  the  alongshore  component  of  momentum. 

The  motion  of  water  inside  the  surf  zone  is  extremely  complex  because  of  the  dy- 
namic upper  and  lower  boundaries  and  generally  three  dimensional  nature  of  the  flow. 
In  order  to  overcome  this  complexity  in  the  surf  zone,  a  number  of  assumptions  have 
to  be  made.  As  a  consequence,  the  models  describing  wave  transformation  and 
longshore  currents  include  empirical  parameters  that  must  be  specified  in  the  formu- 
lations. The  specification  of  these  empirical  parameters  requires  testing  the  model 
against  comprehensive  experimental  data. 

As  ol'  today,  there  is  no  exact  theoretical  formulation  for  the  wave  height  distrib- 
ution inside  the  surf  zone  where  the  waves  break.  Therefore  empirical  formulations  are 
used  to  specify  the  wave  height  in  the  surf  zone.  It  is  generally  assumed  that  the 
breaking  wave  height  in  the  surf  zone  is  a  linear  function  of  the  depth  in  the  form  of: 

U=,h  [1] 


in  which  H  is  the  wave  height,  h  is  the  water  depth  and  y  is  the  proportionality  constant. 
This  relation  holds  fairly  well  for  spilling  breakers.  It  has  been  shown  by  several  inves- 
tigators such  as  Horikawa  and  Kuo  (1967),  Nakamura,  Shirashi  and  Sasaki  (1967), 
Divorky,  Le  mehaute  and  Lin  (1970)  that  this  relation  does  not  hold  for  plunging 
breakers.  Although  this  relation  does  not  hold  for  all  breaker  types  it  has  been  retained 
in  most  models  built  to  predict  the  wave  height  in  the  surf  zone  (Smith  and  Kraus,  1988). 
Wu  et  al.  (1985)  also  utilize  this  relation  in  their  two  dimensional  model  that  includes 
the  bottom  friction  and  the  nonlinear  convective  accelerations.  The  major  deficiency  of 
present  day  models  is  the  inability  to  correctly  describe  wave  breaking  processes. 

B.     OBJECTIVES 

The  objectives  of  this  thesis  are  to  test  and  modify  the  wave  transformation  and 
longshore  current  models  being  considered  for  application  in  the  Tactical  Environmental 
Support  System  (TESS).  These  models  are  part  of  the  the  surf  forecasting  software 
which  was  developed  under  contract  by  Marshall  D.  Earle,  MEC  Systems  Corporation 
to  provide  the  Navy  with  an  improved  Surf  Forecasting  capability.  The  basis  of  this 
model  are  the  wave  transformation  and  longshore  current  models  developed  by 
Thornton  and  Guza  (19S3).  (1986).  In  order  to  make  the  models  more  operational  and 
increase  the  range  o[  applicability,  an  expanded  data  base  of  74  different  data  sets  in- 
cluding field  and  lab  data  are  utilized.  The  model  parameters  are  calculated  by  mini- 
mizing the  error  in  the  least  square  sense  between  the  model  and  the  data,  an  attempt 
is  made  to  predict  the  model  parameters. 


II.     NUMERICAL  MODEL 

The  surf  forecasting  model  predicts  the  wave  height  and  longshore  current  distrib- 
ution across  the  surf  zone.  The  feature  of  this  model  compared  to  most  previous  studies 
is  that  the  waves  are  treated  as  random  rather  than  monochromatic.  The  coordinate 
system  orientation  and  sign  convention  used  in  the  model  is  given  in  Figure  1  .  The 
shoreline  has  been  aligned  with  the  y-coordinate  and  the  x-coordinate  is  positive  off- 
shore. A  wave  orthogonal  approaching  the  beach  from  southerly  direction  has  a  positive 
angle  value  (  +  6)  and  causes  positive  longshore  current,  ie.  to  the  right,  where  as  a 
northerly  wave  orthogonal  is  just  the  opposite. 

A.     DESCRIPTION  OF  WAVE  TRANSFORMATION  MODEL 

The  wave  height  transformation  for  straight  and  parallel  contours  is  calculated  using 
the  energy  flux  balance  in  which  the  onshore  spatial  variation  of  the  energy  flux  is  bal- 
anced by  the  ensemble  averaged  dissipation  due  to  wave  breaking,  <  tb  >  ,  and  the 
frictional  dissipation.  <  tf  >  .   The  energy  flux  balance  is  given  by: 

d(ECgx) 

— t^—  =  <tb  >  +  <lf>  [2] 

where  E  is  the  wave  energy,  C$t  is  the  x-component  of  the  group  velocity. 

Based  on  analysis  of  field  data  acquired  at  Torrey  Pines  Beach,  CA.  it  was  found 
that  the  Rayleigh  distribution  was  in  good  agreement  with  the  data  both  outside  and 
inside  the  surf  zone.  Therefore  the  wave  height  distribution  inside  the  surf  zone  is  de- 
scribed by  the  Rayleigh  probability  density  function  given  by: 

^(//)  =  ^-exr(-(-£-)2)  [3] 

Energy  and  the  x-component  of  the  group  velocity  are  described  using  linear  theory. 
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p{H)dH  =  \pgH2rms  [4] 


y 

i 

i 

^ 

/ 

\                      *•/*. 

f>J» 

2 

\  -e 

.0  ■% 

^' 

-'£** 

/ 

\ 

s  ' 

/ 

\ 

/                    ^  ^ 

B   A 

\ 

/             ^  *** 

E  z 

/     <"  "" 

A   A 

=    / 

\^\ 

H    / 

\         ^\. 

/ 

\                 ^s. 

/ 

\ 

/ 

\ 

/ 

V5     V  O 

/ 

^-^ 

/ 

' 

°%5\ 

A 

r 

*s>s 

~) 

i 

Figure  1.  Orientation  of  the  Coordinate  System  and  Sign  Convention:  The  y-axis 
has  been  aligned  with  the  alongshore  and  the  x-distance  is  positive  off- 
shore.  A  wave  orthogonal  coming  from  the  south  has  a  positive  angle 
value  and  causes  a  longshore  current  to  the  right  as  the  landing  craft 
approaches  from  ofTshore  to  the  beach. 


c    =  <L{1  +      2kh 
°x       2  {        Sinhlkh 


)cos6 


[5] 


where  k  is  the  wave  number  associated  with  peak  frequency/,  h  is  the  water  depth  and 
6  is  the  mean  wave  direction  corresponding  to  the  peak  of  the  wave  spectrum. 

Although  there  are  several  dissipation  mechanisms  (frictional,  percolation,  sediment 
transport,  etc.  )  once  the  waves  start  to  break  turbulent  dissipation  becomes  the  domi- 


nam  dissipative  mechanism.  Therefore  the  model  only  includes  dissipation  due  to  wave 
breaking.  The  frictional  dissipation  is  not  included  in  the  model  because  it  was  found 
that  the  frictional  dissipation  was  less  than  3%  of  the  dissipation  due  to  wave  breaking. 
Even  if  this  is  not  always  true,  the  effect  is  compensated  for  by  the  B  parameter  in  the 
bore  dissipation  mechanism.  It  is  more  reasonable  to  include  the  frictional  effects  in  B 
rather  than  introducing  new  parameters.  Frictional  effects  may  become  significant  in  the 
swash  regions  but  this  is  not  the  concern  of  this  model. 

Unlike  monochoromatic  waves  in  which  there  is  a  single  breaker  line,  the  random 
treatment  of  the  waves  describes  the  larger  waves  to  break  farther  offshore  and  the 
smaller  waves  to  break  closer  to  the  shore  line,  which  is  more  realistic.  Across  the  surf 
zone  at  each  point  there  might  be  broken  and  unbroken  waves  at  the  same  time.  The 
model  calculates  the  breaking  wave  dissipation  utilizing  the  periodic  bore  formulation 
as  suggested  by  LeMehaute(1962).  The  bore  dissipation  mechanism  as  applied  only  to 
the  breaking  waves  is  formulated  as  follows  : 

1  (ft;  -/,,)'  1  IBID3 

E*>«  -  T ""     /,,/,,     Q  *  T « -£~  Q  [«] 

where  p  is  the  water  density,  g  is  the  gravitational  acceleration.  /:,  is  (  h  +  r\  )  ,  /z,  is 
(  h  -  }]  )  ,  ;/  is  the  surface  elevation,  h  is  the  mean  water  depth  and  Q  is  the  volume 
discharge  across  the  bore.  Measurements  done  on  a  gently  sloping  (1:40)  beach  show 
that  the  classical  periodic  bore  formulation  underestimates  the  dissipation  by  30%  to 
(Stive  1983).  The  under-estimation  is  compensated  for  by  adjusting  the  value  of 
the  parameter  B  which  describes  the  amount  of  foam  on  the  crest  of  a  wave  (see 
figure  2).  There  are  several  formulations  suggested  for  the  volume  discharge  across  the 
bore.   The  description  by  Hwang   ei  al.   (1970  )  is  used  in  the  model, 

where  C  is  the  wave  speed.  L  is  the  wave  length  and  h  is  the  mean  water  depth. 

In  order  to  find  the  percentage  of  waves  that  are  breaking,  the  probability  density 
function  is  simply  multiplied  by  a  weighting  funtion. 

Pb(H)  =  W{H)p{H)  [S] 
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Figure  2.  Sketch  of  the  Periodic  Bore  Formulation:  II  is  the  wave  height,  h  is  the 
mean  water  level.  //,  is  the  water  depth  below  wave  trough,  h2  is  the 
surface  elevation  measured  from  the  bottom.  13  is  the  amount  of  foam 
over  the  crest  of  a  breaking  wave. 


Note  that  Pb(U),  is  not  a  pdf  but  a  percentage.  In  order  to  distinguish  the  percentage 
of  the  waves  that  are  breaking,  the  following  empirical  weighing  function  which  fits  the 
observations  is  used  (Thornton  and  Guza,  1983) 
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It  was  found  by  the  same  authors  that  a  value  of  y=.40  and  n  =  2  gave  the  best  fit  to 
the  data. 

The  average  rate  of  energy  dissipation  is  found  by  adding  up  the  dissipation  for  each 
breaking  wave  calculated  using  the  formula 


/        (BID 
ibore  =  -TPg—T— 


[10] 


'I  hen  this  sum  is  divided  by  the  total  number  of  waves  both  broken  and  unbroken.    In 
other  words  the  average  rate  of  energy  dissipation  is  calculated  by  multiplying  the  dissi- 


pation  for  a  single  broken  wave  height  H  by  the  probability  of  wave  breaking  at  each 
height,  as  given  by  Pb{H).   For  the  ensemble 


<zb>=^-Pg{-\     HJPb(II)dH 


[11] 


Substituting  (10)  the  bore  dissipation  is  given  by  : 


[12] 


B.     DESCRIPTION  OF  THE  LONGSHORE  CURRENT  MODEL 

Most  of  the  longshore  current  models  for  straight  and  parallel  contour  beaches  are 
based  on  an  alongshore  momentum  balance.  Assuming  steady  state,  the  alongshore 
momentum  balance  is  given: 


8SV, 
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in  which  the  first  term  on  the  left  hand  side  is  the  wave  induced  momentum  flux,  or  ra- 
diation stress  gradient  and  the  second  term  is  the  lateral  transfer  of  turbulent  momen- 
tum. On  right  hand  side  of  the  equation,  are  the  bottom  shear  and  surface  wind  stresses 
in  the  alongshore  direction. 

The  wave  induced  radiation  stress  due  to  gravity  waves,  first  introduced  by 
Longuet-Higgins  and  Steward  (1964).  is  utilized  to  express  the  change  in  momentum  due 
to  waves.  Thornton  and  Guza  (1983)  using  linear  theory  and  the  energy  flux  equation(2) 
showed  that  the  wave  induced  radiation  stress  gradient  may  be  expressed  as 


55 


[14] 


Neglecting  the  surface  wind  stress  in  the  alongshore  direction  and  substituting  equation 
(12),  the  alongshore  momentum  equation  becomes: 

sin  0„  dS'yx        h 

— cT<t">'— jT-'y  [l5] 

The  tangential  bottom  stress  t*  for  small  angle  of  incidence  and  weak  mean 
longshore  current  can  be  approximated  (  Longuet-Higgins  (1970a,b)  )  : 

t*  =  PC/|u|  F  [16] 

in  which  cf  is  the  bottom  shear  stress  coefficient,  u  is  the  mean  water  particle  velocity 
magnitude  and  V  is  the  mean  longshore  current.  Thornton  and  Guza  (1983)  calculated 
the  mean  water  particle  velocity  magnitude  as  : 


-  Hr„ 


(i)  W 


The  lateral  turbulent  mixing  term  S'yx  is  the  vertically  integrated  Reynold's  stress. 
Since  the  functional  form  of  this  stress  is  not  known,  it  is  usually  parameterized  utilizing 
an  eddy  viscosity  coefficient  v  : 


S'yx  =       pu'v'dz  =  -  pvh  — —  [18] 

The  eddy  viscosity  formulation  by  Longuet-Higgins  (1970b),  which  spatially  decrease 
from  offshore  to  onshore,  is  used  in  the  model  : 

v  =  .Y|.vU£/7  [19] 

where  N  is  an  adjustable  coefficient  and  based  on  dimensional  analysis  has  an  approxi- 
mate range  of  0  <  N  <  0.016    . 

If  the  lateral  transfer  of  momentum  is  not  negligible  and  is  to  be  included,  then  the 
alongshore  momentum  balance  equation  is  a  second  order  ordinary  differential  equation. 
In  this  case,  the  gradient  of  this  term  is  given  by: 


dS'yx 


*lr[M**%]  [20] 


Taking  the  derivatives. 
dS' 


{(**%)  +  (**>*%)+{>*'$■)]   ™ 


Substituting  back  into  the  steady  state  alongshore  momentum  equation,  we  have: 


sin5. 


v  -[,?:  dV   .    ?     7-    dh   dV  ^  ,3/2    </2r 

c>.\  v'£    A      -7-  +  -=-  v 'h  x  —  ——  +  h  '  x - 

d.x        2  dx    dx  rfx2 


-pcf\u\V  [22] 

Grouping  the  terms,  we  end  up  with  the  following  second  order  differential  equation. 


f-v    r   c-r\    1    rI+4"T-* 


rfF 


sinf? 


[23] 


C.     APPLICATION  OF  THE  MODEL 
I.     Wave  Height  Model 

To  apply  the  wave  height  model,  the  the  parameters  B  and  y  are  combined,  re- 
ducing the  degree  of  freedom  to  one.  The  B  parameter  describes  the  amount  of  foam  on 
the  face  of  a  breaking  wave  i.e.  the  intensity  of  wave  breaking,  y  is  a  proportionality 
constant  that  relates  the  rms  breaking  wave  height  to  the  water  depth  and  is  slope  de- 
pendent. The  combined  parameter  is  designed  to  reduce  the  number  of  free  parameters 
to  tune  the  model.  The  term  in  brackets  of  equation  12  in  the  bore  dissipation  equation 
is  a  correction  term  to  the  original  bore  dissipation  formulation.  Therefore  when  the  B 
•  amcters  were  combined  the  y  term  in  the  brackets  was  assigned  a  typical  value 
-  c  imbined  parameter  BG  is  equal  to  —  . 


The  mode]  generated  wave  heights  and  data  were  compared  by  finding  the  BG 
value  that  results  in  the  minimum  root  mean  square  error 


/l     \     {H model"  Hdata)  .--.-. 
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where  Hmode,  is  the  model  generated  wave  height,  Hdcta  is  the  measured  wave  height  and 
n  is  the  number  of  data  points.  The  model  was  iterated  over  the  combined  parameter 
BG  and  the  rms  error  was  computed.  The  BG  parameter  was  varied  from  low  to  high 
values  to  observe  the  model  behavior  and  stability.  Based  on  the  minimum  rms  error, 
the  best  fit  was  obtained.  Examples  from  each  data  set  for  the  model  fit  to  the  wave 
height  data  are  shown  in  Figures  9  to  14. 
2.      Longshore  Current  Model 

The  longshore  current  model  was  compared  with  the  measured  longshore  cur- 
rents for  planar  beaches  and  reasonable  results  were  obtained  with  and  without  turbu- 
lent mixing.  Similar  comparision  for  barred  beaches  showed  that  the  model  fails  to 
predict  the  longshore  current  over  the  shore  side  of  the  bar  when  turbulent  mixing  is  not 
included.  In  order  to  improve  the  results,  turbulent  mixing  parameterized  by  an  eddy 
viscosity  was  included  in  the  model,  which  increased  the  degree  of  longshore  current 
equation  to  the  second. 

This  second  order  differential  equation  (23)  was  solved  numerically  using  a 
centered  difference  scheme.  The  equation  was  put  in  a  finite  difference  form  and  a  linear 
system  of  equations  were  obtained  for  the  interior  of  the  spatial  domain.  At  the 
shoreline,  a  zero  current  boundary  condition  was  assumed  and  the  very  offshore  value 
measured  in  the  data  was  used  as  the  offshore  boundary  condition.  The  system  of 
equations  was  put  in  a  tridiagonal  matrix.  The  tridiagonal  matrix  was  solved  by  a  con- 
venient tridiagonal  solver  by  Maron  (1987).  A  mid  value  of  0.008  was  assigned  to  the 
lateral  mixing  coefficient,  N  ,  which  has  an  approximate  range  of  0.0  to  0.016  and  the 
model  was  iterated  over  the  bed  shear  coefficient  until  the  error  between  the  model 
generated  current  and  the  data  is  minimized.  Then  with  the  best  shear  stress  coefficient 
obtained  the  model  was  iterated  over  N  from  0.0  to  .030  until  the  rms  error  was  reduced. 
Examples  of  model  fits  of  Superduck  and  Santa  Barbara  data  are  shown  in  Figure  15  and 
16  respectively. 


III.     DATA  AND  .METHODS 

A.     WAVE  HEIGHT  TRANSFORMATION 

To  test  and  modify  the  wave  transformation  model,  74  data  sets  were  used.  The 
data  consisted  of  both  laboratory  and  field  data.  These  data  provided  many  different 
wave  conditions  from  the  long  period,  low  steepness  Pacific  swell  to  steeper  Atlantic  and 
North  Sea  storm  waves.  The  experiments  were  conducted  over  various  bottom  profiles 
including  both  planar  and  barred  beaches.  The  data  sets  are  described  below.  The  data 
sets  are  named  after  either  the  site  of  the  data  or  the  person  that  acquired  the  data.  In 
the  tables  that  describe  the  data,  tan/?  is  the  mean  beach  slope,  f  is  the  peak  frequency 
in  hertz.  H.„:  is  the  root  mean  square  wave  height,  —r-  is  the  deep  water  wave  steepness, 
c0  is  the  deep  water  surf  similarity  parameter,  BG  is  the  combined  model  parameter  and 
e..s  is  the  rms  error  percentage  between  the  model  generated  wave  height  and  the  data. 
The  deep  water  surf  similarity  parameter 


-^  [27] 


was  calculated  for  each  data  set  in  which  //.  is  the  deep  water  wave  height,  Lc  is  the  deep 
.    .  and  tan/;  is  the  mean  beach  slope  (Battjes  1974).    The  deep  water 

wave  height  was  calculated  from  the  measured  wave  height  in  shallower  water  using  the 
linear  wave  theory  shoaling  coefficient  corresponding  to  the  the  peak  frequency.  Simi- 
larly, the  deep  water  wave  length  was  calculated  from  the  peak  frequency  using  linear 
wave  theory. 

1.     Torrey  Pines  (  NSTS  )Data 

Torrey  Pines  wave  height  data  were  acquired  as  part  of  the  Nearshore  Sediment 

Transp  -  STS I  at  Torrey  Pines  Beach.  CA.  during  November  1978.  The  beach 

'..mar  and  homogeneous  in  alongshore  direction.    Shadowing 

effects  of  the  offshore  islands  and  refraction  limits  the  angle  of  wave  incidence  in  the  surf 

zone  to  a  narrow  band  of  less  than  +  15  degrees  and  usually  much  less.    The  winds 

z  light  and  variable  during  the  experiment.  An  extensive  array 


of  instruments  was  deployed  to  study  the  nearshore  wave  dynamics.  The  sensors  used 
in  this  study  were  11  Marsh-McBirney  electromagnetic  current  meters,  four  Stathem 
pressure  transducers  and  four  resistance  wire  wave  staff.  The  wave  heights  were  deter- 
mined by  the  zero-up  crossing  method,  in  which  the  wave  height  is  defined  as  the  dif- 
ference between  maximum  and  minimum  within  two  consecutive  zero-up  crossings.  The 
data  were  filtered  to  include  only  the  wind  wave  band  of  frequencies  between  0.05  and 
0.3  Hz  .  The  data  are  described  in  Table  1.  Details  of  the  experiment  are  given  in 
Thornton  and  Guza  (1983). 

Table  1.     TORREY  PINES  WAVE  HEIGHT  DATA 


CASE 

DATA 

TYPE 

tan/? 

f 

(Hz) 

Hfms 

(cm.) 

HQI 1^ 

L 

BG 

(%) 

Nov  04 

Field 

.021 

.063 

35.0 

.0007 

.78 

11.8 

7.1 

Nov  10 

Field 

.021 

.055 

56.0 

.0008 

.72 

13.6 

10.9 

Nov  12 

Field 

.021 

.077 

81.0 

.0027 

.39 

7.6 

4.0 

Nov  17 

Field 

.023 

.069 

47.1 

.0011 

.69 

13.8 

7.3 

Nov  18 

Field 

.023 

.069 

53.0 

.0012 

.65 

11.2 

3.8 

Nov  20 

Field 

.023 

.063 

56.2 

.0011 

.71 

12.2 

5.9 

2.  Santa  Barbara  (NSTS)  Data 

The  Santa  Barbara  data  (Table  2)  were  acquired  during  the  second  NSTS  ex- 
periment conducted  at  Leadbetter  Beach  from  30  January  to  23  February  1980.  The 
site  was  selected  because  the  waves  are  homogeneous  in  the  alongshore  direction  and  the 
bottom  contours  are  relatively  straight  and  parallel.  The  Leadbetter  Beach  is  composed 
of  wellsorted  fine  to  medium  size  sand.  The  mean  nearshore  slope  varied  between  .03 
to  .06  during  the  experiment.  The  shoreline  has  an  unusual  east-west  orientation  along 
a  predominantly  north-south  coast.  The  open  ocean  waves  are  limited  to  a  narrow 
window  of  approach  -9  degrees  to  +9  degrees  centered  at  249  degrees.  Therefore  the 
waves  from  the  open  ocean  were  directionally  narrow  banded  swell  arriving  at  relatively 
large  angles  of  about  10  degrees  to  the  shoreline.  Instrumentation  included  24 
Marsh-McBirney  electromagnetic  current  meters  and  14  pressure  sensors.  The  exper- 
iment is  described  in  detail  in  Thornton  and  Guza  (1986). 

3.  Battjes  and  Stive  (1985)  Data 

The  first  three  cases  (  Table  3  )  are  lab  data  that  were  acquired  in  a  wave  flume 
using  mechanically  generated  random  waves.   The  bottom  profiles  include  a  schematized 


Table  2.     SANTA  BARBARA  WAVE  HEIGHT  DATA 


CASE 

DATA 

TYPE 

tan/) 

f 

(Hz) 

(cm.) 

HJL, 

Z. 

BG 

erms 
(%) 

Eeb  02 

Field 

.059 

.063 

32.4 

.0010 

1.8 

5.6 

10.2 

Feb  03 

Field 

.1  '-- 

.070 

54.6 

.0023 

.q 

s,  j 

12.0 

Per  (H 

Field 

.038 

.070 

56.0 

.0024 

7 

6.6 

".8 

Feb  05 

Field 

.035 

.078 

45.1 

.0023 

.7 

4.S 

7.5 

Feb  06 

Field 

.033 

.090 

26.4 

.OOP 

.8 

8.4 

8.8 

bar  profile  and  a  plane  slope  in  concrete  (Case  5  and  Case  10).  Case  15  has  a  profile 
with  bar  and  trough  on  sand.  All  of  the  lab  wave  height  data  were  measured  by  means 
of  parallel  wire  conductivity  wave  gages. 

Case  IS  data  is  field  data  acquired  on  a  beach  near  Egmond  City.  The 
Netherlands.  The  field  site  was  chosen  to  have  more  or  less  statistically  uniform  condi- 
tions alongshore  and  normal  wave  incidence.  These  field  data  were  acquired  under 
stormy  conditions  during  which  the  incident  rms  wave  height  was  2.8  meters.  It  was 
noted  that  the  wind  influence  was  negligible.  The  Hrmj  was  calculated  from  the  surface 
elevation  time  series  assuming  that  the  surface  elevation  is  a  narrow  band  Gaussian 
process,  which  appears  to  be  in  good  agreement  with  observations  (  Thornton  and  Guza 
1983  i.  The  rms  wave  height  is  related  to  the  standard  deviation,  o.  of  the  surface  ele- 
vation bv 
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Table  X     BATT.JES  AND  STIVE  WAVE  HEIGHT  DATA 


CASE        ^j™|    tan/? 

f 

(Hz) 

Hrmi 
(cm.  1 

//   L 

- 

BG 

(  %  ) 

i     l-y-     '     .«»34     '     .443 

11.8 

.016 

.2"      1      3.1 

8.7 

1     Lab        .025    |    .663     |     13.6    | 

.13 

2." 

15  " 

15        i     Lab 

.012 

.557 

13.2     1       ."29 

.07 

1.0 

7.5 

i    Field 

.010 

.115 

2"^     1       .026 

.07 

2.2 

12.6 

4.     Superduck  (1986)  Data 

The  Superduck  data  (Table  4)  were  acquired  during  the  one  month  long  ex- 
periment held  by  the  U.S.  Army  Corps  of  Engineers  at  the  Field  Research  Facility  in 
Duck,  N.C.  in  October  19S6.  The  experiment  was  called  Superduck.  The  site  of  the 
experiment  is  along  an  extensive  barrier  island  formation  known  as  North  Carolina's 
Outer  Banks.  There  are  no  major  structures  along  the  shoreline.  The  beach  has  a  mean 
foreshore  slope  of  1:10.  usually  a  single  bar  formation  and  a  mean  slope  of  1:100  off- 
shore the  bar.  Extensive  bathymetric  data  were  acquired  using  the  Coastal  Research 
Amphibious  Buggy  (CRAB)  so  that  the  wave  height  and  longshore  currents  could  be 
processed  with  accurate  and  detailed  bathymetry  data.  The  waves  and  currents  were 
measured  using  a  pressure  sensor  and  a  vertical  array  of  three  Marsh- McBirney  current 
meters  mounted  on  a  moveable  sled.  The  details  of  the  data  can  be  found  in  Whitford 
and  Thornton  (1989). 


Table  4. 

SUPERDUCK  WAVE  HEIGHT  DATA 

CASE 

DATA 
TYPE 

tan/? 

f 
(Hz) 

H,ms 
(cm.) 

#o/£o 

£0 

BG 

(%) 

Oct  15 

Field 

.037 

.160 

92.0 

.<>207 

25 

3.4 

9.7 

Oct   16 

Field 

.049 

.  1 90 

9S.1 

.0239 

.26 

2.S 

13.4 

Oct  17 

Field 

.041 

.190 

72.0 

.0175 

.31 

9.4 

8.1 

Oct  IS 

Field 

.038 

.190 

93.0 

.0226 

2S 

5.2 

8.2 

5.     Vincent  (1984)  Data 

The  Vincent  data  (Table  5)  were  acquired  from  a  series  of  experiments  done  at 
the  Coastal  Engineering  Research  Center  (CERC)  in  a  wave  tank  of  44  meters  long.  0.6 
meters  wide  and  a  water  depth  of  0.6  meters.  At  one  end  of  the  tank  was  a  smooth 
concrete  ramp  with  a  slope  of  1:30  .  A  random  wave  generator  was  used  to  generate 
waves.  The  wave  height  data  were  measured  using  nine  resistance  wave  gauges  (Vincent, 
19S4). 


Table  5. 

VINCENT  WAVE  HEIG 

HT  DATA 

CASE 

DATA 
TYPE 

tan/? 

r 

(Hz) 

11,,., 
(cm.) 

HJU 

*. 

BG 

e,mi 
(%) 

1058 

Lab 

1:30 

.885 

8.6 

.0456 

16 

.6 

10.8 

1116 

Lab 

1:30 

.763 

11.2 

.0452 

16 

.8 

5.9 

1132 

Lab 

1:30 

.769 

13.1 

.05  3  S 

14 

.8 

7.6 

1148 

Lab 

1:30 

.671 

14.6 

.0461 

17 

.8 

5.0 

1206 

Lab 

1:30 

.671 

15.7 

.0499 

16 

.8 

6.0 

1509 

Lab 

1:30 

.403 

6.2 

.0065 

41 

2.6 

2.9 

1526 

Lab 

1:3<» 

.403 

9.3 

.0098 

34 

2.2 

6.6 

1604 

Lab 

1:30 

.403 

12.1 

.0128 

30 

1.6 

6.9 

1623 

Lab 

1:30 

.403 

14.5 

.0153 

27 

1.8 

5.6 

1638 

Lab 

1:30 

.403 

15.2 

.0160 

27 

1.6 

7.5 
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Figure  ?.  Setup  In  C  LRC  Wave  Tank:  Ten  gages  were  placed  three  of  them 
were  located  over  the  horizontal  bottom  and  seven  of  them  over  the  1:30 
slope. 


6.     CERC  (1983)  Data 

Wave  height  data  (Tables  6  and  7)  were  acquired  in  Coastal  Engineering  Re- 
search Center  (CERC)'s  wave  tank.  The  tank  was  0.9  meters  high,  0.4  meters  wide  and 
45.7  meters  long.  The  setup  of  this  tank  was  designed  to  obtain  data  on  the  shoaling 
and  breaking  random  waves.  Ten  parallel  resistance  wire  wave  gages  were  used.  Three 
of  them  were  located  over  the  horizontal  bottom  of  the  tank  and  seven  of  them  over  the 
1:30  concrete  slope  (  Figure  3  ). 


Table  6. 

CERC  WAVE  HEIGHT  DATA  (  PART  I  ) 

CASE 

DATA 
TYPE 

tan/? 

f 
(Hz) 

Hrm, 

(cm.) 

H0ILq 

tB 

BG 

(%) 

30 

Lab 

1:30 

.699 

6.1 

.0199 

24 

4.1 

6.7 

31 

Lab 

1:30 

.758 

6.3 

.0241 

21 

4.1 

6.7 

32 

Lab 

1:30 

.667 

6.7 

.0204 

23 

2.7 

2.2 

33 

Lab 

1:30 

.505 

3.6 

.0064 

42 

3.3 

5.7 

34 

Lab 

1:30 

.757 

6.8 

.0261 

21 

2.1 

1.9 

35 

Lab 

1:30 

.6S4 

3.8 

.0121 

30 

3.5 

5.0 

36 

Lab 

1:30 

.684 

7.6 

.0242 

21 

2  2 

1.9 

37 

Lab 

1:30 

.Soo 

6.4 

.0270 

20 

2.4 

2.4 

38 

Lab 

1:30 

.709 

6.5 

.0259 

21 

2.3 

2.1 

39 

Lab 

1:30 

.684 

4.4 

.0139 

28 

3.6 

5.1 

4() 

Lab 

1:30 

.657 

6.7 

.0200 

24 

2.7 

3.0 

41 

Lab 

1:30 

.564 

9.0 

.02ol 

24 

2.1 

2.0 

42 

Lab 

1:30 

.558 

4.4 

.0096 

34 

3.9 

4.6 

43 

Lab 

1:30 

.284 

4.7 

.0023 

69 

5.0 

3.5 

44 

Lab 

1:30 

.289 

9.8 

.0051 

47 

2.6 

6.9 

45 

Lab 

1:30 

.598 

6.7 

.0167 

2o 

2.6 

3.4 

40 

Lab 

1:30 

.320 

7.4 

.0049 

48 

3.6 

6.0 

47 

Lab 

1:30 

.316 

5.8 

.0037 

54 

4.7 

5.1 

4S 

Lab 

1:30 

.452 

5.7 

.0081 

37 

3.9 

4.3 

49 

Lab 

1:30 

.680 

6.5 

.00:4 

->-} 

2.2 

2.1 

50 

Lab 

1:30 

.6SO 

6.5 

.O2o4 

23 

2.1 

2.1 

51 

Lab 

1:30 

.598 

6.7 

.0167 

26 

2.7 

3.4 

Table  7. 

CERC  WAVE  HEIGHT  DATA  (  PART  II  ) 

CASE 

DATA 
TYPE 

tan/? 

f 

(Hz) 

(cm.) 

M)/^0 

L 

BG 

errrs 
(%) 

52 

Lab 

1:30 

.452 

5.S 

.0082 

37 

3.9 

4.3 

53 

Lab 

1:30 

.680 

6.5 

.0204 

23 

2.2 

2.1 

54 

Lab 

1:30 

.680 

6.5 

.0204 

23 

2.1 

2.2 

55 

Lab 

1:30 

.598 

6.7 

.0167 

26 

2.7 

3.6 

56 

Lab 

1:30 

.543 

7.1 

.0147 

27 

3.3 

1.8 

57 

Lab 

1:30 

.239 

7.7 

.0025 

66 

4.7 

4.1 

58 

Lab 

1:30 

.234 

7.7 

.0024 

68 

4.0 

5.6 

59 

Lab 

1:30 

.4b  1 

5.7 

.0085 

36 

4.4 

3.0 

60 

Lab 

1:30 

.370 

5.7 

.0052 

46 

4.8 

1.5 

61 

Lab 

1:30 

.384 

3.8 

.0037 

54 

7.0 

2.8 

62 

Lab 

1:30 

.384 

3.8 

.0037 

54 

7.0 

2.6 

63 

Lab 

1:30 

.543 

7  2 

.0148 

27 

3.4 

2.0 

6*1 

Lab 

1:30 

.680 

9.8 

.030S 

19 

2.0 

3.4 

65 

Lab 

1:30 

.680 

9.8 

.0308 

19 

1.9 

3.3 

66 

Lab 

1:30 

.598 

9.8 

.0243 

21 

2.4 

1.9 

67 

Lab 

1:30 

.5C?S 

9.9 

.0247 

21 

2.3 

1.5 

6S 

Lab 

1:30 

.281 

8.6 

.0042 

52 

3.3 

4.7 

6  1 

Lab 

1:30 

.277 

8.6 

.0040 

53 

3.3 

3.6 

70 

Lab 

1:30 

.294 

10.2 

.0056 

45 

3.2 

3.9 

"1 

■  Lab 

1:30 

.649 

5.4 

.0157 

27 

3.7 

3.8 

72 

Lab 

1:30 

.667 

6.1 

.0187 

24 

3.3 

2.8 

73 

Lab 

1:30 

.598 

5.6 

.0140 

28 

3.9 

3.0 

"4 

Lab     |     1:30 

.602 

5.9 

.0148 

27 

3.7 

3.0 

B.     LONGSHORE  CURRENTS 

The  longshore  current  model  described  in  Chapter  2  is  applied  to  the  Santa  Barbara 
and  Superduck  data  sets.  These  are  the  only  high  quality  longshore  current  data  avail- 
able acquired  on  a  homogeneous  alongshore  bathymetry  with  adequate  measure  of  di- 
rectional wave  forcing.  The  Santa  Barbara  beach  profiles  were  almost  planar  while  the 
Superduck  profiles  were  barred.  Of  the  combined  two  months  of  data  collection,  five 
cases  from  Santa  Barbara  and  four  cases  from  Superduck  were  selected  for  comparision. 

Since  the  longshore  currents  are  induced  by  the  wind  waves,  simultaneous  meas- 
urements of  the  wave  field  and  longshore  current  measurements  is  essential.  This  re- 
quirement was  satisfied  in  the  data  sets.  In  the  Tables  8  and  9  describing  the  longshore 
currents,  cf  is  the  bed  shear  stress  coefficient,  N  is  the  lateral  mixing  coefficient,  f  is  the 
frequency  in  hertz.  Hrms  is  the  the  rms  wave  height,  —r2-  is  the  deep  water  wave  steepness, 
t0  is  the  deep  water  surf  similarity  parameter,  erm  is  the  rms  error  percentage  between 
the  model  generated  longshore  current  and  the  data. 


Table  8.     SANTA  BARBARA  LONGSHORE  CURRENT  DATA 


CASE            cf 

N 

f 
(Hz) 

(cm.) 

HJLo 

e 

L 

(%) 

Feb  02 

.0085 

.0160 

.063 

32.4 

.0010 

6.4 

1.8 

50.0 

Feb  03 

.0090 

.0190 

.070 

54.6 

.0023 

7.8 

.90 

5  5 . 5 

Feb  04 

.OU90 

.0001 

.070 

56.0 

.0024 

9.0 

.77 

27.9 

Feb  05 

.0000 

.0010 

.078 

45.1 

.0023 

8.4 

.73 

22.1 

Feb  06 

.00  SO 

.0020 

.090 

26.4 

.0017 

8.3 

.79 

31.0 

Table  9.     SUPERDUCK  LONGSHORE  CURRENT  DATA 

CASE 

cf 

N 

f 

(Hz) 

H0 

(cm.) 

HJL, 

e 

L 

ermi 
(%) 

Oct  15 

.0025 

.OolO 

.160 

92.0 

.0207 

9.1 

.25 

5.9 

Oct  16 

.0035 

.0030 

.  1 90 

98.1 

.0239 

14.7 

.26 

4.5 

Oct  17 

.0065 

.0001 

.190 

72.0 

.0175 

6.3 

.31 

7.1 

Oct  18 

.0005 

.0070 

.  1 90 

93.0 

.0220 

11.8 

.25 

17.7 

IV.     DISCUSSION 

The  wave  transformation  model  was  applied  to  74  different  data  sets.  The  very 
offshore  wave  height  and  direction  measured  were  used  as  an  input  to  the  model.  The 
model  is  not  very  sensitive  to  the  exact  orientation  of  waves  to  the  shoreline  since  the 
group  speed  computations  involves  the  cosine  of  the  angle  of  wave  incidence  and  the 
cosine  of  a  small  angle  is  not  significantly  different  from  one.  Therefore  exact  wave 
orientation  is  not  required  and  the  error  due  to  wave  orientation  is  expected  to  be  small. 

The  data  sets  and  the  profiles  were  compiled  and  the  model  was  iterated  over  the 
BG  parameter  to  reduce  the  root  mean  square  error  percentages  in  the  least  square 
sense.  A  fully  developed  bore  would  be  expected  to  have  a  B  value  of  1  and  since  the 
typical  proportionality  constant  y  for  Rayleigh  distribution  is  about  0.4.  the  combined 
parameter  BG  for  a  fully  developped  bore  would  be  equal  to  6.25.  The  statistics  of  the 
BG  parameter  for  the  whole  74  data  sets  were  computed.  A  mean  of  BG  parameter  of 
4  with  a  standard  deviation  of  2.9  was  obtained  for  all  the  data  sets. 

Model  sensitivity  to  the  choice  of  the  y  parameter  was  investigated  by  Gill  (19S5). 
He  varied  the  y  parameter  from  0.35  to  0.55  and  found  that  the  error  increase  was  less 
than  8-10%.  In  the  same  manner,  the  B  parameter  was  iterated  from  0.9  to  1.7  and  it 
was  found  that  the  increase  in  rms  error  was  less  than  10-12%.  Similar  type  of  investi- 
gations were  done  by  Thornton  and  Guza  (1983)  for  Torrey  Pines  data  set  and  it  was 
found  that  varying  the  B  parameter  +  25%  resulted  in  an  increase  of  rms  error  less  than 
10%. 

The  combined  parameter  BG  is  plotted  against  the  rms  error  in  Figure  5,  and  it  is 
seen  that  changing  the  BG  parameter  one  standard  deviation  either  side  of  the  mean 
causes  an  increase  in  rms  error  of  less  than  9%  .  Based  on  these  investigations  it  is 
evident  that  the  model  is  not  overly  sensitive  to  the  exact  choice  of  parameters. 

Relatively  high  values  of  BG  parameter  are  required  to  minimize  the  rms  error  for 
Torrey  Pines  data.  The  Santa  Barbara  data  had  almost  the  same  deep  water  steepness 
but  didn't  require  as  high  values  of  BG.  The  difference  is  hypothesized  due  to  the  defi- 
ciencies of  linear  bore  theory,  for  spilling  type  breaker  the  turbulence  is  confined  to  a 
surface  layer,  basically  between  the  crest  and  trough  of  a  wave.  For  plunging  or  surging 
type  breaking  the  turbulence  is  not  confined  between  the  crest  and  trough,  but  pene- 
trates much  deeper.    The  actual  depth  h,  concerning  the  volume  discharge  computation 
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Figure  4.      Variation  of    B    with  Deep  Water  Surf  Similarity  Parameter:      The  B 

parameter  shows  a  linear  trend  with  the  deep  water  surf  similarity  pa- 
rameter, t0.  The  dashed  line  is  the  linear  regression  line.  Notice  that  the 
B  parameter  values  increase  as  the  breaker  type  shifts  from  spilling  to 
plunging  {t0  >  0.4). 


will  be  different  from  the  mean  water  level  as  expressed  by  the  linear  bore  theory.  It  is 
evident  that  there  is  deviation  from  the  bore  formulation.  Although  the  incident  wave 
parameters  such  as  deep  water  steepness  and  frequency  for  Torrey  Pines  and  Santa 
Barbara  data  sets  are  about  the  same,  the  beach  slope  of  Santa  Barbara  is  about  twice 
that  of  Torrey  Pines.  This  result  implies  that  the  wave  breaking  is  a  function  of  both 
deep  water  wave  steepness  (~-  )  and  beach  slope  (tan/?).  Therefore  the  wave  dissipation 
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Figure  5.  Variation  In  Rms  Error  Duo  to  Change  In  BG  Parameter:  Changing 
BG  (BETAGAMMA)  parameter  one  standard  deviation  on  either  side 
of  the  mean  does  not  increase  the  error  greater  than  9%. 


rate  is  lower  over  mild  slopes  compared  to  steep  slopes.  'I  he  combination  of  plunging 
type  breaker  and  the  beach  slope  effect  is  thought  to  be  responsible  for  high  values  of 
BG  for  Torrey  Pines  data  set.  This  idea  is  confirmed  when  the  BG  parameter  obtained 
for  all  cases  plotted  versus  deep  water  steepness.    As  seen  in  (Figure  6),  as  the  deep 


water  steepness  becomes  less  than  0.02  the  BG  parameter  tends  to  increase  exponen- 
tially. For  deep  water  steepness  greater  than  0.02  the  BG  parameter  stays  invariant  at 
a  value  of  about  2.9.  The  curve  fitted  to  the  individual  points  has  the  form  of 
y  =  c.  +  —  where  o.  is  equal  to  2.25  and  /?  is  equal  to  0.0085. 

Using  the  deep  water  surf  similarity  parameter  the  breaker  type  is  classified  as  fol- 
lows (Battjes,  1974)  : 

£fl  <  0.4  Spilling 

0.4  <  l0  <  2.0  Plunging 

C„  >  2.0  Surging  or  Collapsing 

Most  of  the  data  sets  tabulated  in  chapter  2  fall  in  the  spilling  range  with  a  few  of  them 
being  in  the  plunging  range.  There  were  no  surging  or  collapsing  type  breaker.  Torrey 
Pines  and  Santa  Barbara  (NSTS)  data  differ  from  the  other  data  due  to  their  low  steep- 
ness, about  one  order  of  magnitude  smaller  than  the  other  data  sets.  As  a  result,  the 
breaker  type  based  on  the  surf  similarity  parameter  by  Battjes  (1974)  suggest  plunging 
for  all  cases  of  the  NSTS  data  sets.  Breaking  intensity  increases  from  spilling  to  plunging 
to  collapsing  breakers  as  <fc  increases.  Therefore  it  is  expected  B  and  hence  BG  will  be 
a  function  of  t0  and  increase  as  c0  increases. 

The  BG  parameters  obtained  for  all  data  sets  are  plotted  versus  deep  water  surf  pa- 
rameter and  an  exponential  trend  was  observed.  A  curve  in  the  form  of y  =  ae"*  is  fitted. 
Excluding  one  outlying  data  point  seen  in  the  right  side  of  Figure  7  ,  it  is  found  that 
a  =  1.303  and  /?=  2.4  fitted  the  data  well.  This  curve  acts  in  a  very  similar  manner  to  the 
breaking  wave  height  curves  in  the  Shore  Protection  Manuel,  U.S.  Army  Corps  of  En- 
gineers (1984)  in  this  exponentially  increasing  region  of  the  curve.  Utilizing  this  func- 
tion, the  model  parameter  can  be  predicted.  The  deep  water  surf  parameter  can  be 
computed  for  a  particular  beach  of  interest  and  inputting  this  value  to  the  fitted  curve 
function,  the  BG  parameter  can  be  predicted. 

The  BG  parameter  was  plotted  versus  the  incident  rms  wave  height.  No  correlation 
was  evident.  Therefore  it  was  concluded  that  the  model  parameter  is  not  strongly  de- 
pendent on  incident  wave  height  by  itself.  Running  all  the  data  sets,  it  was  observed  that 
the  model  is  capable  of  handling  a  small  scale  wave  height  of  6  cm  lab  data  as  well  as  a 
27S  cm  larse  wave  heieht  measured  during  a  stormv  weather. 
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Tigure  6.  DC;  Versus  Deep  Water  Steepness:  The  solid  curve  represents  the 
curve  fitted  to  BG  values  versus  deep  water  steepness  for  all  data  sets  in 
the  form  of  j  =  a.  +  —. 
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Figure  7.  BG  Versus  Deep  Water  Surf  Parameter:  The  solid  curve  represents 
the  curve  fitted  to  BG  values  versus  deep  water  surf  parameter  for  all 
data  sets  in  the  form  of  y  =  o>.ep*. 
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Figure  8.  Rms  \Va\e  Height  and  Depth  Profile  of  Sample  Torrey  Pines 
Data:  Model  predicted  wave  height  versus  November  20  1978  data. 
The  solid  line  is  the  model  predicted  wave  height  while  the  individual 
points  are  the  measured  wave  heights.    1  he  second  plot  shows  the  depth 

profile. 
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Figure  9.  Kins  Wave  Height  and  Depth  Profile  of  Sample  Santa  Barbara 
Data:  Model  predicted  wave  height  versus  February  03,  1980  data. 
The  solid  line  is  the  model  predicted  wave  height  while  the  individual 
points  are  the  measured  wave  heights.  1  he  second  plot  shows  the  actual 
depth  profile. 


Sallenger  ei  al.  (19S5)  investigated  the  dependence  of  wave  height  on  local  slope  and 
depth.   They  found  that  y  linearly  increases  with  bottom  slope  in  the  form  of: 

V  =  3.2  tan  fi  +  0.30  [29] 

They  also  concluded  that  the  y  coefficient  was  independent  of  the  deep  water  wave 
steepness.  The  value  of  y  was  computed  for  all  the  data  sets  using  the  above  relation. 
The  coefficients  obtained  based  on  this  relation  ranged  from  0.33  to  0.48  with  a  mean 
of  0.40  .  Their  conclusion  that  y  was  independent  of  the  deep  water  wave  steepness  was 
found  also  to  be  valid  for  the  present  data  sets  tabulated  in  chapter  3.  The  B  values  were 
computed  by  substituting  the  calculated  y  values  into  the  combined  BG  parameter.  The 
computed  B  parameter  values  in  this  manner  ranged  from  0.46  to  1.24  with  a  mean  of 
0.S2  and  a  standard  deviation  of  0.17.  This  result  shows  that  B  is  assigned  reasonable 
values  and  agrees  with  the  physics  of  the  problem.  These  B  parameter  values  showed 
an  increasing  trend  with  the  deep  water  surf  similarity  parameter  (see  Figure  4)  .  As  the 
breaker  type  shifts  from  spilling  to  plunging  (;,,  >  0.4),  the  B  parameter  values  increase. 
The  longshore  current  model  was  run  first  neglecting  turbulent  mixing.  Relatively 
high  rms  error  percentages  in  the  least  square  sense  were  obtained  (See  Figure  15  ). 
No  matter  how  much  c,  was  adjusted,  it  was  not  possible  to  obtain  longshore  current 
over  :he  shore  side  of  the  bar  where  the  waves  reform  and  are  no  more  breaking.  Based 
on  this  result,  it  was  decided  to  include  turbulent  mixing  in  the  model.  VlcDougal  et 
al.  1 1986)  investigated  the  influence  of  lateral  mixing  on  longshore  currents.  They  tested 
a  variety  of  lateral  mixing  models  and  concluded  that  the  longshore  current  profile  is 
insensitive  to  the  form  of  the  mixing  model.  Therefore  it  was  decided  to  use  the  lateral 
mixing  model  of  Longuet-Higgins  (1970b)  which  has  the  form  of  v  =  Nxjgh  .  Inclusion 
of  the  turbulent  mixing  increased  the  degree  of  the  longshore  current  equation  to  the 
second.  More  reasonable  results  are  obtained  by  including  the  turbulent  mixing  term 
(Figure   15). 
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Figure   10. 


Rms  Wave  Height  and  Depth  Profile  of  Sample  Datljes  and  Stive 
Data:  Model  predicted  wave  height  versus  Case  IS  data.  The  solid 
line  is  the  model  predicted  wave  height  while  the  individual  points  are 
the  measured  wave  heights.  The  second  plot  shows  the  actual  depth 
profile. 
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Figure   11. 


Rms  Wave  Height  and  Depth  Profile  of  Sample  Supeiduck 
Data:  Model  predicted  wave  height  versus  October  16,  19S6  data. 
The  solid  line  is  the  model  predicted  wave  height  while  the  individual 
points  are  the  measured  wave  heights.  The  second  plot  shows  the 
depth  profile. 
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Figure  12.  Rms  Wave  Height  and  Depth  Profile  of  Sample  Vincent  Data:  Model 
predicted  wave  height  versus  Case  1509  data.  The  solid  line  is  the 
model  predicted  wave  height  while  the  individual  points  are  the  meas- 
ured wave  heights.  The  second  plot  shows  the  constant  slope  bottom 
profile. 
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Figure  13.  Rms  Wave  Height  and  Depth  Profile  of  Sample  CERC  Data  Part 
(I):  Model  predicted  wave  height  versus  Case  44  data.  The  solid  line 
is  the  model  predicted  wave  height  while  the  individual  points  are  the 
measured  wave  heights.   The  second  plot  shows  the  depth  profile. 
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Figure  14.  Rms  Wave  Height  and  Depth  Profile  of  Sample  CERC  Data  Part 
(II):  Model  predicted  wave  height  versus  Case  68  data.  The  solid  line 
is  the  model  predicted  wave  height  while  the  individual  points  are  the 
measured  wave  heights.   The  second  plot  shows  the  depth  profile. 
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Figure   15.      Longshore  Current  and  Depth  Profile  of  Sample  Superduck  Data:      1  he 

dotted  line  indicates  the  model  generated  longshore  current  without 
turbulent  mixing.  The  solid  line  is  the  model  predicted  current  with 
turbulent  mixing  and  the  individual  points  arc  the  measured  wave 
heights.  The  second  plot  shows  the  depth  profile.  Notice  that  the  ve- 
locities on  shore  side  of  the  bar  significantly  change  depending  on 
inclusion  or  exemption  of  the  turbulent  mixing. 
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Longshore  Current  and  Depth  Profile  of  Sample  Santa  Barbara 
Data:  I  he  dotted  line  indicates  the  model  generated  longshore  current 
without  turbulent  mixing.  The  solid  line  is  the  model  predicted  current 
with  turbulent  mixing  and  the  individual  points  are  the  measured  wave 
heights.    The  second  plot  shows  the  depth  profile. 


In  the  Santa  Barbara  data,  the  model  generates  a  maximum  longshore  current  closer 
to  the  shoreline  compared  with  the  maximum  current  measured  (see  Figure  16).  This 
result  was  regardless  of  the  inclusion  or  exemption  of  the  turbulent  mixing.  Therefore, 
the  inclusion  of  turbulent  mixing  to  the  model  did  not  cause  significiant  improvement 


for  Santa  Barbara  longshore  current  data.  The  model  comparison  was  found  insensitive 
to  the  N  value.  Based  on  this  result,  turbulent  mixing  may  be  excluded  in  the  model  for 
planar  beaches.  As  disscussed  earlier  this  is  not  true  for  barred  beaches.  This  would 
reduce  the  number  of  model  parameters  to  be  specified  and  also  since  the  longshore 
current  equation  will  stay  in  the  first  order,  some  cpu  time  will  be  saved.  Planar  beaches 
are  more  appropriate  for  Amphibious  Warfare  compared  to  barred  beaches.  Landings 
will  not  be  made  on  barred  beaches  unless  necessary.  If  a  barred  beach  is  to  be  chosen 
for  landing,  then  the  turbulent  mixing  should  be  included. 

The  mean  value  of  the  bed  shear  stress  coefficient.  cf  for  the  present  longshore  cur- 
rent data  is  0.006  and  the  mean  lateral  mixing  coefficient  is  0.0055.  Since  there  is  not 
enough  longshore  current  data  to  obtain  predictive  relations  for  the  model  parameters, 
these  mean  values  are  suggested  for  use  in  model  implementation.  Longshore  current 
magnitude  is  directly  proportional  to  cr  It  should  be  recognized  that,  since  the  optimal 
cf  values  varied  between  0.0025-0.009,  the  expected  current  magnitude  can  have 
significiant  error  using  the  suggested  values.  For  generality,  since  it  does  not  matter  to 
include  or  exclude  the  turbulent  mixing  term  for  planar  beaches,  the  model  can  run  with 
these  c,  and  N  values  regardless  of  the  beach  profile.  These  values  should  be  checked 
as  more  longshore  current  data  becomes  available. 


V.     SUMMARY  AND  CONCLUSIONS 

In  this  study,  the  wave  height  transformation  and  longshore  current  models  of 
Thornton  and  Guza  (1983,  1986)  are  applied  to  a  variety  of  wave  conditions  over  planar 
and  barred  beaches  for  the  purposes  of  testing  and  evaluation.  74  data  sets  of  wave 
height  transformation  and  9  data  sets  of  longshore  currents  are  utilized.  The  model  re- 
sults are  optimized  in  the  least  square  sense  to  obtain  the  best  fit  to  the  data.  The  model 
coefficients  in  the  wave  transformation  model,  y  a  measure  of  saturated  wave  breaking 
and  B  a  measure  of  breaking  wave  intensity,  are  combined  into  a  single  BG  parameter. 
Based  on  data  analysis,  the  combined  BG  parameter  in  the  wave  transformation  model 
is  found  to  be  correlated  with  the  deep  water  surf  parameter.  Therefore,  this  correlation 
can  be  used  to  predict  the  single  model  BG  parameter  in  future  model  applications. 

Testing  the  model  using  field  and  laboratory  data  shows  that  the  wave  transforma- 
tion model  is  highly  robust  in  describing  the  wave  height  distribution  inside  the  surf 
zone.  Varying  the  combined  BG  parameter  one  standard  deviation  from  the  mean  of  the 
74  data  sets,  did  not  increase  the  rms  error  greater  than  9%.  Therefore  it  is  concluded 
that  model  is  not  overly  sensitive  to  the  parameters. 

The  y  parameter  was  calculated  using  the  relationship  by  Sallenger  and  Holman 
(1985)  and  the  B  values  were  calculated  from  the  BG  parameter.  The  mean  propor- 
tionality coefficient  y  for  the  whole  74  data  sets  was  0.4.  The  B  parameter,  which  de- 
scribes the  intensity  of  wave  breaking,  appears  also  to  be  correlated  with  the  deep  water 
surf  parameter.  This  is  expected  since  the  deep  water  surf  parameter  is  used  to  predict 
breaker  type,  and  B  is  a  measure  of  breaking  wave  intensity. 

lor  longshore  current  computations  the  turbulent  mixing  needs  to  be  included  over 
barred  beaches,  but  for  planar  beaches  its  inclusion  has  an  insignificant  effect  on  the 
results.  Therefore  it  is  suggested  the  turbulent  mixing  be  included  in  a  general  model  for 
consistency.  At  present  there  is  not  a  means  for  predicting  the  appropriate  model  pa- 
rameters for  longshore  current  computations.  Therefore  it  is  suggested  the  longshore 
current  model  be  run  with  the  mean  values  obtained  for  the  bed  shear  stress  and  lateral 
mixing  coefficients,  0.006  and  0.0055  respectively,  regardless  of  the  beach  profile. 

The  model  is  operationally  handy  since  it  is  only  initialized  by  the  offshore  wave 
height,  direction  and  period,  which  can  be  observed  aboard  the  ship.    The  construction 


of  the  model  is  quite  simple,  and  does  not  require  much  computation  time.    Therefore 
the  model  can  run  on  a  personal  computer  with  almost  instant  response. 
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